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Abstract 



Two interaction mechanisms of particles in a fluid are proposed on base 
of forces, mediated by hydrodynamic thermal fluctuations. The first one is 
similar to the conventional van der Waals interaction, but instead of been 
mediated by electromagnetic fluctuations, it is mediated by fluctuations of 
hydrodynamic sound waves. The second one is due to a thermal drift of 
particles to the region with a bigger effective mass, which is formed by the in- 
volved surrounding fluid and depends on an inter-particle distance. The both 
mechanisms likely can be relevant in interpretation of the observed long-range 
attraction of colloidal particles, since a set of different experiments shows the 
attraction energy of the order of fcgT and, perhaps, only a fluctuation mech- 
anism of attraction provides this universality. 

PACS numbers: 82.70.Dd, 05.40.+j 
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I. INTRODUCTION 



Charge-stabilized colloidal suspensions exhibit a variety of unusual physical properties 
[ Q-0]. Colloidal particles can be organized into crystal M and into structures with clus- 
ters and voids p|-|T0|j. A system of colloidal particles may undergo different types of phase 
transitions fTT|-|TT||. The topological phase transitions in two dimensional systems of col- 
loidal particles have been reported in Refs. |IEj,I5|. In Ref. PO] buckling instabilities in a 
confined colloidal crystals were analyzed. An interesting behavior of colloids in external 
fields were reported in Ref. |^TJ. Colloidal particles accept in an electrolyte some effective 
charge screened by counterions at Debye's length \e>, what is described by the repulsion 
potential of Derjaguin, Landau, Verwey, and Overbeek (DLVO) The DLVO theory, 

as a result of solution of the linearized Poisson-Boltzmann equation, has been questioned 
in Refs. [|l2, p2^ . The generalization of DLVO interaction via a modification of counterion 
screening was reported in Ref. JJ3] . 



Despite a long history of the problem, the interaction of colloidal particles remains a 



matter of challenge and controversy. The authors of Ref. [|12|| , studying phase transitions 
in charged colloidal systems, found a substantial deviation from predictions resulting from 



a screened Coulomb interaction. Experimental data [21] suggest a net attraction of par- 
ticles for explanation of measurements (see also the Comment on that work P^fl). Later 
an attraction of micron size particles separated by a micron-sized distance has been estab- 
lished experimentally. The interaction potential has been found to have a minimum —Uq at 
center-to-center particle distance Ro of the micron size. In the work [2£§ colloidal particles, 
confined between two glass planes, exhibit Uq ~ 0.2&gT for different ionic strengths. This 
type of experimental arrangement, studied in another lab [p7|1 , gives Uq ~ (0.3 — O.^fc^T 
for different particles diameter and distances between planes. Another study of colloidal 
particles, confined between two glass planes, gives Uq ~ 1.3/cbT or less for different par- 
ticle diameters and distances between planes [0,^9[ (see also [pO| , pTl ). The attraction of 
0.5k B T per neighbor particles (separated by the micron-sized distance) has been deduced 
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for the colloidal crystal in Ref. |30||. Polystyrene colloidal particles of a diameter 0.5/xm on 



the water-air interface exhibit Uq ~ 0.5ksT at R Q ~ 0.9/mi [p2| . For micron sized colloidal 
"molecules" on the air- water interface the binding energy can be estimated as a few kT. 



For bound particles on the fluid-fluid interface an approximate estimation is Uo ~ AkT [|33 
As one can see, despite of different conditions (even particles on interfaces) there is a very 
stable common feature of all various experiments: the attraction minimum is always of the 
order of ksT. This provides a hypothesis of some common mechanism of the micron-sized 
attraction, which is responsible for the universality of Uq ~ &bT. Note, mean-field energies 
in colloidal physics (electrostatic and hydrodynamic) are a few orders of magnitude bigger 
comparing to kT at room temperature. 

Let us analyze some mechanisms of attraction proposed in the literature. 

A principal question is that can like-charged particles separated by a micron distance 
attract each other due to a solely electrostatic mean-field interaction (for example, by some 
charge redistribution) which is not accounted by DLVO theory. This type of attraction 
was predicted in Ref. |33| . The results of Ref. are not applicable to dielectric particles 
but only to those containing inside the electrolyte identical to outside one. The correct 
calculation of the interaction of that kind shows only repulsion in contrast to Ref. . The 



electrostatic attraction between like-charged particles based on the mean-field approach 
was also predicted in Ref. f35fl , but that conclusion was incorrect, as shown in the works 



36| , |37[| . An attraction potential force of like-charged particles at a micron distance due to a 



mean-field mechanism seems to be extremely unlikely and the works [0,0] provide strong 
arguments in support of this statement. 

The conventional van der Waals attraction, mediated by the high-frequency (visible light) 
electromagnetic fluctuations, is very small in the micron range, consisting less than 10 _2 T 



at room temperature p8|-^l|1 , and cannot provide the observed attraction. Measurements of 
this type of attraction in colloidal systems at short distances of the order of a few hundred 
Angstroms have been performed in Refs. H2|fl3]. An analysis of attraction at such short 



distances, including influence of the surface roughness, is given in Refs. [|4J-f47 



The van der Waals interaction mediated by the low frequency (order of the plasma 
frequency) electromagnetic fluctuations decays fast with a distance R between two particles 
as — Vo exp(— 2R/Xd) PH^TOf . Vq has in our case the thermal fluctuation origin and is less 
than strong electrostatic repulsion energy even at R ~ Ad. The exponential dependence 
exp(— R/Xd) of the attraction potential in Ref. dealing with the same effect, differs 



from the above correct exponent [f|8],£|0]]. The results {|9[ cannot explain the attraction at 



R ~ IOXd in Ref. [^8] at least for the reason the attraction, reported in [f49 |, is too small on 
that distance. 



The effect of correlation of counterions [p0|-|5^] (see also [pq -p7||) results in their re- 
arrangement in the vicinity of a colloidal particle of the order of mean distance between 
counterions, which is less than Xr>. On that short distance an attraction is possible as 



shown in Refs. 0-54|. From a position on the long micron-sized distance, the correlation 



of counterions forms the effective charge Ze of the particle, which is screened at Ad and 
determines the repulsive part of the interaction. 

Depletion forces between colloidal particles originate from an influence of the thermody- 
namic energy by the finite size of small particles, which constitute the surrounding medium 
The range of the depletion interaction is determined by the size of small parti- 
cles, which can be not necessary "small" like in Ref. |yj, where surrounding polymer coils 
have the macro-particle size. When the surrounding electrolyte contains only microions, the 
depletion force is of a very short range comparing to the micron scale. 



In Ref. [§3| a mechanical effect is proposed for explanation of the observed motion of 
particles towards each other, when they move away a single wall due to the Coulomb repul- 
sion |30| , |31| . This mechanical effect is irrelevant in analysis of interaction of free particles in 
a fluid or ones confined between walls or on interfaces (see also ||64|| ). 

As one can conclude, no one above mechanism is responsible for the set of observations 



of long distance attraction of colloidal particles (may be excepting the mechanical effect p5 



which can be relevant for the particle motion in the one wall geometry j3U|). An important 
hint in search and selection of possible mechanisms of attraction is the universality of a 



depth of the attraction minimum Uq ~ ksT. A most probable mechanism, satisfying this 
criterion, is an interaction, mediated by thermal fluctuations of some physical values. This 
is a microscopic type of interaction, which can be formulated in terms of a potential energy. 
Regardless of specificity of fluctuating matter, the free energy of thermal fluctuations is 
always proportional to ksT. For example, the attraction potential of two dielectric particles 
(refractive index n) in the water (refractive index no) due to thermal fluctuations of the 
electromagnetic field is proportional to c(n,n )kBT. If n is close to no, then c ~ {n — no) 2 , 
according to the perturbation theory. With real values of n and no for polystyrene and water 
the coefficient c is very small, what makes the conventional van der Waals attraction in the 
micron range to be negligible. Nevertheless, if to put formally n — > oo (zero fields gradient on 



the particles surface), the attraction in the micron range is not small comparing to ksT [38 
This is similar to the situation in hydrodynamics, where fluctuating electromagnetic waves 
are substituted by hydrodynamic ones and the fluid velocity is zero on the particles surface. 
Hence, one can expect the interaction, mediated by fluctuating hydrodynamic waves, to be 
not small like in the electromagnetic case with n — > oo. 

The nature of forces, mediated by fluctuations of sound waves in a fluid, can be clarified 
in the following way: the energy of thermal fluctuations of the fluid depends on a distance 
between two particles, which play a role of obstacles for a fluid motion, and hence this results 
in a force. The non-electromagnetic type of fluctuation forces, discussed by Dzyaloshinskii, 
Lifshitz, and Pitaevskii | 65[] , can be considered as some sort of the van der Waals []38| 



or 



the Casimir |66| forces. We do not concern here the story of origin of the two titles for 



fluctuation interactions. 

In addition to the attraction, mediated by fluid fluctuations, there is another mechanism 
of a fluctuation interaction related to hydrodynamics: this is an interaction, mediated by 
thermal fluctuations of particle positions in a fluid. The effective particle masses depend on 
a fluid mass, involved into the motion. The fluid mass depends on an inter-particle distance 
and therefore the effective particle masses also depend on that distance. If a mass of a classic 
non-dissipative harmonic oscillator depends on the coordinate, the mean coordinate shifts 



to the region with a bigger mass (smaller velocity), since the particle spends there more 
time. Analogously, for a Brownian motion of particles in a fluid a thermal drift occurs in 
the direction of a bigger effective mass. As shown in the paper, the effective particle mass, 
which determines the effect, should be calculated on base of the Euler (non-dissipative) 
hydrodynamics and it can be called the Euler mass. 

The above fluctuation interactions relate to hard spheres (no mean-field interaction ex- 
cepting an infinite repulsion on contact). If particles are charged, the fluctuation hydro- 
dynamic forces are weakly modified by the Coulomb effects in a bulk fluid under typical 
experimental conditions. The only interaction, mediated by fluctuations of particle posi- 
tions, can be influenced by an electrostatic restriction of particle fluctuation motions in the 
direction perpendicular to charged confining plates. The total interaction is a sum of the 
fluctuation ones and the Coulomb repulsion. 

The following results are presented in the paper: (i) the attraction forces between two 
parallel plates and spherical particles with short inter-surface distance, mediated by hydro- 
dynamic fluctuations of sound waves, are shown to exist and they are calculated analytically; 
(ii) a novel interaction of particles in a fluid is proposed, which is based on dependence of 
their effective masses on a distance between them; (iii) the Fokker-Planck equation for two 
particles in a fluid is derived; (iv) a difference in measurements of a particle interaction by 
optical tweezers and by a long time statistics is pointed out and calculated. 

II. THE HYDRODYNAMIC VAN DER WAALS INTERACTION 

Suppose two particles are placed inside a hydrodynamic medium, they are totally fixed 
in space, and serves only as obstacles for a fluid motion. There is no macroscopic motion 
in the system and the only motion is caused by thermal fluctuations of the fluid velocity 
v(f, t). In this case the free energy of thermal fluctuations of the fluid F(R) depends on the 
distance R between bodies. The function 




(1) 
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is an interaction mediated by fluid fluctuations. Analogously to the conventional van der 
Waals interaction mediated by electromagnetic fluctuations, the potential (|1|) can be called 
the hydrodynamic van der Waals interaction. To find the free energy of thermal fluctuations 
of the fluid one can start with the linearized Navier-Stokes equation |57[ 

p Yt = ~ Vp + vV2g + + 3 ) vdivt7 ^ 

There are two types of fluid motion, one of them is a transverse diffusion and the second 
one is longitudinal sound waves, associated with the density variation. The equilibrium free 
energy of transverse motions is determined by the Boltzmann distribution of their kinetic 
energies and does not depend on the friction coefficient in the thermal limit. Since there is no 
static interaction for transverse motions, their equilibrium free energy in the thermal limit 
depends on the total volume, but not on relative positions of bodies. Therefore, transverse 
fluctuations do not result in an interaction. Quite opposite situation occurs for longitudinal 
motions, when the total free energy is a sum of energies of different sound modes. The 
spectrum of sound waves depends on the distance between bodies R due to hydrodynamic 
boundary conditions on body surfaces and this results in independence of the free energy. 
Hence, the fluctuation interaction between bodies is mediated by hydrodynamic sound waves 
like the conventional van der Waals interaction is mediated by fluctuations of electromagnetic 
ones. Putting v = Vd<p/dt, one can obtain from Eq. (g) 

Through thermodynamic relations and the continuity equation one can obtain 5p = 
— PSqV 2 0, where so is the adiabatic sound velocity |)7|. At the typical frequency uj ~ s /a 
(a is the particle radius), involved into the problem, the dissipative term in Eq. (|j) is small 
and one can write 

<H - agV 2 * = (4) 



dt 2 

According to the case of a small friction, the boundary condition V n = to Eq. 
corresponds to the Euler equation |)7] . From general point of view, free energy of a system 



of harmonic oscillators does not depend on friction in the thermal limit. 

Let us consider first the case of two infinite parallel plates, separated by the distance R, 
when the frequency spectrum has the form ||67|| . 

i 2 n 2 
R 2 



u\(k) = sUk 2 + ^-) (5) 



The free energy per unit area of the system now can be expressed as a sum of energies of 
independent oscillators 

^/^f>^> (6) 

J ^ > n=l 



This type of interaction has been considered in literature |39| . The simplest way to calculate 
the energy (|6|) is to divide the whole interval R by small segments ao = R/N. Then 

n=l v / n=1 n=1 \ / 

Using the relation a;]^[^° =1 (l + x 2 /7c 2 n 2 ) = sinhx, the Stirling formula for N\, omitting 
constants, and i?-linear terms, one can obtain from Eqs. (|§-[7p 

Now one should put in Eq. (Bp Na^ = R and after integration we obtain, according to Eq. 
(H) , the van der Waals interaction per unit area of two infinite plates 

= "If tp < 9) 

Despite the van der Waals interaction (H) formally coincides with the electromagnetic van 
der Waals formula at e — » oo [Q, existence of such effect in hydrodynamics is not trivial as 
follows from Section III. 

The result (^j) enables to calculate the van der Waals energy U v dw(R) °f t wo spheres 
of radii a, separated by the center-to-center distance R, when (R — 2a) a. In this case 
the inter-surface distance, measured in the direction parallel to the center-to-center line, is 



S 



R — 2a + (x 2 + y 2 )/a, where x and y are the coordinates in the plane perpendicular to the 
center-to-center line. The procedure turns out to an integration 

f ( x 2 + y 2 \ 
U vd w{R) = dxdy u vdW [R - 2a H — J (10) 

The result of integration at (R — 2a) <C 2a is 

Uvdw{R) = -^T-^— (spheres) (11) 

lo R — la 

The interaction energy (|ll|) is calculated under fixed boundaries of particles, which cor- 
responds to zero sound velocity s of the particles material (infinite acoustic mismatch). 
Under reduction of the mismatch U v dw decreases, turning to zero at s — Sq, if there is 
only longitudinal acoustic mode inside particles. Sound velocities for polystyrene particles 
s ~ 2.1 x 10 5 cm/s and for water s — 1.5 x 10 5 cm/s provide a finite mismatch, which reduces 
the result (|TTD, as one can show, approximately four times. But in reality this conclusion is 
not correct, since transverse acoustic modes of the particle material increase the mismatch 



turning the interaction towards the result (|TT|). An account of the finite mismatch is a 



matter of further study, but the exact U v dw seems to be close to the result (|TT|). 

The above calculations are applicable, strongly speaking, to uncharged particles like hard 
spheres. For charged particles an electric charge density of a fluid en is finite and decays 
on the Debye length A/> It results in a modification of the spectrum of fluctuations by 
adding the plasma frequency s 2 k 2 + uj 2 n/no, where u 2 = Airne 2 /eM (M is the mass of the 
fluid molecule) and uq ~ 10 23 cm -3 is the molecular density of the fluid. The coefficient 
n/riQ comes from that the sound spectrum is formed by the all fluid, but the Coulomb 
gap is determined only by a small amount of ions n ~ Z/AiraX 2 ^ ~ 10 17 cm~ 3 (a typical 
Z ~ 10 4 ). Since the wave vector k is inversely proportional to the inter-particle distance, 
one can conclude from here that the Coulomb effects modify the above hard sphere result 
(p!l|), when (R — 2a) exceeds 100/mi. 
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III. THE PARADOX 



In case of the conventional van der Waals interaction, mediated by electromagnetic fluc- 
tuations, the mean values of electric and magnetic fields are zero < E >=< H >= 0. The 
stress tensor cr^ for electromagnetic field is quadratic with respect to fields and hence the 
mean value < cr^ > is not zero. This makes the origin of the force due to electromagnetic 
fluctuations to be straightforward. The situation with hydrodynamic fluctuation forces is 
different. For the linearized Navier-Stokes equation the mean value of the velocity is zero 
< v >= (the same for a fluctuation part of p). The hydrodynamic stress tensor 

°ik = v (lr + Sr) ( 12 ) 



\dr k dri 

is linear in fluctuation variables, its fluctuation part is zero, and a fluctuation force has to be 
zero in this approximation. A non-zero contribution to < > can result from the nonlinear 
terms in the Navier-Stokes equation neglected in the above approach. Account of this non- 
linearity has been done in Ref. |7(J and a finite fluctuation force has been obtained. This 



result was shown to be incorrect in Ref. JTJj , where the exact mean value of the stress tensor 
(|T2|) was found to be zero on the base of exact non-linearity of the Navier-Stokes equations. 
A conclusion of Ref. is that a hydrodynamic fluctuation interaction is impossible, which 
contrasts with the result ([|). What is on a wrong track? 

To understand the situation let us consider the linear chain of small particles, connected 
by elastic springs, shown in Fig. [I] and described by the dynamic equation 

d 2 u n s 1 



dt 2 b 2 



(u n+ i + m„_i - 2u n ) (13) 



where s is the sound velocity and b is the period. Two no-moved big particles substitute 
small particles, as shown in Fig. [l|. The system is elastic and the force acting on a big 
particle, placed on the site n is 

I \ /1A\ 

F n = -p2-{u n +l ~ Un-l) (14) 

°0 
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Here m is the mass of a a small particle. The free energy of fluctuation motion of 
small particles U V( jw is determined by summation on self-frequencies uJi of the system 
F = \.n(Hui/T), resulting in 

U vdW = ^lnW (15) 

where N is a number of springs between two big particles. Two different positions of big 
particles, "natural" in Fig. |l](a) and "compressed" in Fig. 0(b), have identical self frequencies 
since the system is harmonic and hence U^ w = uf} w (equal N in Eq. (|15|)). In this situation 
there is no van der Waals force, what is clear, since the mean value of the linear force ([T4|) 
should be zero. On the other hand, for another "natural" position in Fig. |I](c) self frequencies 
differ from those in Figs. |l](a) and 1(b), in this case U^ w ^ U^ C J W (different N in Eq. (|T5|)) 
and the van der Waals force is non-zero. It is a consequence of that a transition from the 
position (a) to the position (c) in Fig. [I] cannot occur within a harmonic approximation, 
one should destroy some harmonic springs and rearrange them again in a different way. The 
linear expression for force ( |H| ) is not valid to describe a transition from (a) to (c) and a real 
force is non-linear, what makes its average finite even for < u n >= 0. 

An analogous situation takes place in hydrodynamics. According to its derivation, Eq. ([|) 
is valid only for discrete R = Nbo, where bo is inter-atomic distance and N is an integer 
number. The full dependence of u V( iw on R has a structure on the atomic scale corresponding 
to remove of discrete atomic layers from the inter-plane region. The hydrodynamic expres- 
sion for the stress tensor (|l|) is not valid at such short scale, like Eq. (|14D cannot describe 
breaking of harmonic bonds. In contrast to smooth van der Waals potential, mediated by 
electromagnetic fluctuations, the interaction, mediated by hydrodynamic ones, has a struc- 
ture as a function of distance on the atomic scale, superimposed on the smooth function ([|). 
To some extend, this is analogous to observation of the structured interaction potential |62 



where the role of atoms was played by small particles. The resulting statement is that the 
hydrodynamic expression for the stress tensor (TT2|) cannot be used for calculation of fluctu- 
ation forces since it becomes non-linear (and contributes to those forces) at short distances 
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where the hydrodynamic approach is not valid. Hydrodynamic fluctuation forces should be 

on 



calculated on the base of energy like it is done in this paper. The conclusion of Ref. |7T 
absence of hydrodynamic fluctuation forces based on use of the hydrodynamic stress tensor 
is incorrect. 

IV. INTERACTION MEDIATED BY FLUCTUATIONS OF PARTICLE 

POSITION 

The hydrodynamic van der Waals interaction is formed by fast density fluctuations with 
a typical frequency ojl ~ s /a ~ 10 9 s _1 . Besides this longitudinal motion, there are also 
slow transverse fluctuations of the fluid driven by fluctuations of particle linear velocities u\ 
and «2 of the small frequency ujt ~ r]/pa 2 ~ 10 6 s _1 . This type of transverse fluctuations 
is independent of longitudinal one. Thermal fluctuations of particle velocities mediate an 
important part of the total interaction, which supplements U V( iw- 

Before consideration of this contribution we focus at first on some aspects of derivation of 
non-linear dissipative equations, which are useful to understand a formation of interaction 
by fluctuations of particle velocities. Let us start with a one-dimensional motion in the 
potential V(x) of a particle with the mass m(x), attached to a heat bath, which provides a 
friction. A convenient way to proceed is to use the formalism of Caldeira and Leggett |72 
of the infinite set of thermally equilibrium oscillators. The Lagrangian has the form 



mix) .0 
— — x 



where F(x) = J x dz^r](z) is a non-linear coupling to the thermostat. Using the formalism of 
Caldeira and Leggett, one can derive the Langevin equation in the limit of high temperatures 

1 dm i 

m{x)x + -- i 2 + r](x)x + V'(x) = y/r](x)f(t) (17) 

where the average is defined as (f(t)f(t')) = 2T5(t — t'). Let us suppose the viscosity r\ 
sufficiently big, which separates the big frequency Q ~ rj/m from the low frequency V"/r) of 
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the viscous motion in the potential V(x). One can also derive from the Lagrangian ([16]) the 
Fokker-Planck equation for the distribution function W, which in the low frequency limit 
depends only on x and t 

dW(x,t) d f W d r- , . r/ T dW] 



dt dx \ rj{x) dx rj{x) dx 

The additional potential in Eq. fll8|) is 



T 

I(x) = — — - lnm(x) (19) 



The Langevin ( |TTD and the Fokker-Planck ( |18D equations are derived independently from 
the initial system ( p]| ) and do not contain uncertainties since the high frequency limit is well 
defined by Eq. ([H]). If to omit two mass terms in Eq. QT7|), the high frequency limit becomes 
indefinite (Q = oo) and an attempt to derive the Fokker-Planck equation from the Langevin 



equation meets the Ito-Stratonovich uncertainty [73J as a result of loose of a definition of 
the high frequency limit. Note, neither the Ito nor the Stratonovich approaches result in the 
correct Fokker-Planck equation (|18D for a massive particle with a variable viscosity rj(x). 

To clarify an origin of the effective potential (|1^) let us represent the variable x in Eq. (|T7D 
as x(t) + Sx(t), where the small correction Sx(t) varies rapidly with frequencies Q ~ r]/m and 
x(t) is a slow varying variable. In Eq. ( |TTD one can keep the second order of 8x, considering 
x(t) as an instant argument. After an average over high frequencies Eq. ([17]) turns over into 
a low frequency part with two fluctuation induced terms 



T)(x)x + 



3x 2 dx {6X) 2^{x) 3x {6Xf) 



Vv@)f(t) (20) 



By means of the fluctuation-dissipation theorem [7|] one can write (see also |75] 



71 J -oo muJ + lr l m 



The integration path in Eq. ( ]2T| ) can be deformed into the far semi-circle at the upper 
half-plane of the complex u. Hence, (Sx 2 ) is determined by high frequencies and the use of 
instant variables, leading to incorporation of (5x 2 ) into the low frequency equation (^0|), is 
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correct. Analogously 



, fc/)= _M P r_^_ = i (22) 

' ' _ J-oo ^{muj + irj) y/fj 



7T 



The letter P means the principal value of the integral, which can be represented as the inte- 
gration along the infinite contour, consisting of the real axis plus a small circle around zero 
at the upper half-plane, and minus the integration along that small circle. The integration 
along the infinite contour can be shifted to the far upper half-plane and gives zero, but the 
integration along the small circle gives a finite result. Hence, (Sxf) is determined by zero 
frequency limit, the use of instant argument and incorporation of (Sxf) into Eq. (|20D are not 
appropriate. Among two fluctuation induced term in Eq. (|20|) only one (with dm/dx) has 
sense and it results in the potential ([19]). This remains true in a more general case, when 
the linearized form of the dynamic equation is 

(-mco 2 - iur)(uj) + V")5x w = (23) 

In this case the mass m is strictly defined by the condition r)(u)/u — > at uj — > oo and the 
potential ( |i~9"D is determined by the high frequency limit according to Eq. (|2~ID , when only a 
mass term plays a role. 

An origin of the potential (0) can also be also understood from the following non- 
rigorous arguments. One can write formally the free energy F = —T\n(ApAx/h), where 
momentum fluctuations (Ap) 2 ~ mT are formed on a short time scale Q~ l and fluctuations 
of the coordinate (Ax) 2 ~ Tt/rj are slow. In the expression 

T T h 2 n 

F = lnm+-ln— J- (24) 

2 2 TH v ' 



the first term originates from fast fluctuations of momentum and corresponds to Eq. (|19|). 

The interaction (|T9"| ) has a simple interpretation. Suppose the classical non-dissipative 
particle with the variable mass m(x) moves with the total energy E in the harmonic potential 
ax 2 . For a variable mass the mean displacement (x) ^ 0, since in the region with bigger 
mass the particle spends more time having a smaller velocity. One can easily show, when 
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m(x) varies slowly on the scale of the particle amplitude, (x) can be calculated putting the 
constant mass m(0) and adding the potential — (E/2) lnm(x). In other words, a particle 
tends to be more time in a region with bigger mass. This conclusion remains correct for a 
dissipative case with fast fluctuations of velocity instead of harmonic oscillations, when the 
energy E should be approximately substituted by the temperature T. This correspond to 
the potential fll9D . 

One can easily generalize the method for multi-dimensional case. Suppose in the high 
frequency limit the kinetic energy has the form 

K = ^ rriijifyRiRj (cu -> oo) (25) 

Then in the equation of motion, where only the kinetic part is kept, 

mjUVi, + - ~ ^) *A = Fi (26) 

one can consider the variables again as a sum of slow and fast parts Ri (t) +5Ri (t) . An average 
of the quadratic (with respect to 5R) part produces the fluctuation induced effective force 

F * f = F + l^ (SRkSRj) (27) 
With account of the average (5RkSRj) = Tm^f the effective force reads 



2 dRi ' dR t 
From here the generalization of Eq. fllTf) follows 



T 

— ln(detm) 



(28) 



I(R) = ~ln detm(.R) (29) 

The additional potential energies (p!9|) or (p9[) have a fluctuation origin and are mediated 
by fast fluctuations of velocity with the typical frequency Q ~ rj/m, which form, for each 
instant coordinate, the quasi equilibrium free energy. One can formulate a rule of calculation 
the effective fluctuation potential even for a system with a complicated dynamics: one has 
to find the kinetic energy (p5]) in the limit of high frequency and to insert the mass tensor 
into Eq. (|2"5|). The method in the presented form is applicable to classical systems with the 
kinetic energy proportional to square of velocities. 
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V. TWO PARTICLES IN A FLUID AND THE EULER MASS 



Suppose there are two particles in a fluid of radii a separated by the center-to-center 
distance R. If they perform an oscillatory motion with a high frequency u, the fluid velocity 
obeys the Euler equation everywhere in a fluid excepting a thing layer of the thickness 



~ y^rj/pcu close to the particle surfaces, where the full Navier-Stokes equation should be 
used |J7J. Hence for finding the mass tensor ( fffi) one has to solve the Euler equation with 
the boundary condition for a normal component of the fluid velocity. For this reason, the 
mass, corresponding to the high frequency limit of particle dynamics, can be called the Euler 
mass. For example, the Euler mass tensor of one particle is 

rriij = —a (2p + p)Oij (30) 

where p is the mass density of the particle. See also Refs. f)8| ,|69||. 

In the case of two particles in a bulk fluid the Euler mass cannot be calculated analytically 
for an arbitrary relation between R and a. Nevertheless, there is a situation, when an 
analytical calculation of I{R) in a full range is possible. This is a case of two particle 
confined between two parallel plates separated by the distance of the particle diameter, 
which corresponds to the current experiments ||80|| . Let us suppose the particles to be of 
cylindrical shape with the axis of the length 2a perpendicular to the plates. In this case 
particle and fluid velocities are directed parallel to the plates and the problem becomes two 
dimensional. A velocity of the incompressible Euler fluid can be written as v = V0 with 
a boundary condition for the normal derivative dcj)/dn canceling the normal components 
of particle velocities. Since the scalar cf) satisfies the Laplas equation V 2 = in two 
dimensions, one can use a conformal transformation to convert the geometry into a planar 
one. If the centers of two particles localize at positions Kez = ±R/2, Im^ = at the complex 
z-plane, the conformal transformation 

R 2 ,\ 1/2 . w 



-a 2 



coth — (31) 
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maps the two circles on two infinite parallel lines 



, i R R 

Re w = ± In h \ — - 

2a V 4a 2 



(32) 



at the complex w-plane. For the plane geometry the problem can be solved in a straightfor- 
ward way. The mass tensor is easily diagonalyzed by two center-of-mass motions and two 
relative ones. Using the formula (^9j), one can obtain after some calculations 

2(R 2 - 4a 2 ) ^ n(R - VR 2 - 4a 4 ) 2 



/(it!) 



Tin 



-Tin 



1 - 



E 



\2n 



1 + Po/p) ^i(R+ VR 2 - 4a 2 ) 2 ™ + (2c 



\2n 



1 + 



2(R 2 - 4a 2 



E 



n 



(R - VR 2 - 4a 2 ) 



2/i 



a 2 (l + po/p) ^ (i2 + V^R 2 - 4a 2 ) 2 ™ - (2c 
At po = p in limiting cases one obtains 



\2n 



(33) 



I(R) 



ln(7r 4 /72) + (6/ir 2 )y/(R-2a)/a ; {R - 2a) < 2a 



(34) 



-a 4 //? 4 ; 2a < i? 

J(i?) is plotted in Fig. |^ by a dashed line. The van der Waals interaction for cylindrical 
particles of the length 2a can be obtained from Eq. (|9|) in the same manner like Eq. (|10|). 
At (R — 2a) -C 2a for two cylinders of the length 2a 



U v . 



dW 



C(3) 



T 



3/2 



(cylinders) 



(35) 



16 \R-2a / 

The total interaction U v dw{R) + I{R) is shown in Fig. |2| by the solid line, where the close 
particles limit (|35l) is extrapolated up to the region (R/2a — 1) < 1. The interaction I(R) is 
not influenced by the Coulomb effects in a bulk fluid, since it is mediated by fluctuations, 
corresponding to an incompressible fluid. Nevertheless, in a restricted geometry, for example 
between two charged plates, fluctuation motion in some direction can be restricted, which 
can modify the interaction I(R). 



VI. THE FOKKER-PLANCK EQUATION FOR TWO PARTICLES IN A FLUID 



The low frequency dynamic equations of particles 



QiRH.+cnRy.^F, 



?■>■ 

1,2 



(36) 
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is non-linear due to a coordinate dependence of friction coefficients in Eq. (|36|). The (j-term 
in Eq. |36| results in a macroscopic force on the particle (1) due to a macroscopic motion of 
the particle (2). This is called the hydrodynamic interaction [[T6|-|79|. 

The Fokker-Planck equation for the distribution function W(R,t) of two particles in a 
fluid can be derived in a manner like one leading to Eq. (|T8D. The result is the following 

dW d 



Of 



-i 



dU tot (R) BW 
W — + T 

■ dR 



OR, 



j j 



(37) 



where Qj = ((q — (l J )/2 and the total interaction potential consists of some mean- field 
electrostatic part U (R) and two fluctuation interactions 



U to t(R) = U{R) + U vdW (R) + I(R) 



(38) 



The friction tensor in the limit a -C R has the form |76 

3a 9a 2 \ 



3na 



+ 



IR ' 1QR 2 )" v ' \AR ' 16i? 2 J R 2 
For colloids experiments, which deal with long time statistics [^,^], the equilibrium dis- 
tribution function W ~ exp(—Utot(R) /T) is relevant. In the optical tweezers experiments 
|28H3T] two particles are initially fixed and statistics of initial motions after release is stud- 
ied. If friction coefficients would be coordinate independent, the both methods give the 
same interaction potential. But situation becomes different for coordinate dependent fric- 
tion coefficients. If at the moment t — the distribution function was artificially localized 
by optical tweezers at the point Rq (W = 5(R — Rq)) the average inter-particle distance 
(Ri) = j d 3 RiW after release at t = obeys the relation 





Su + 



27a 2 ^ RiRj 



(39) 



dt 



(Ri) = Q^Ff 



(40) 



where the effective force is 



dUtot Tr d x 
' dRi Up 8R q ^ 



(41) 



This force consists of the potential part and the noise-induced drift, it is measured in the op- 
tical tweezers instant experiments in contrast to the long-time-statistics experiments, which 



give only the first potential term. Under experimental conditions |29| the noise-induced drift 
is small. But if two particles are electrostatically fixed in the middle plane between two glass 
plates, separated by the length h (R <C h), the effect of the noise-induced drift on the initial 
motion is not small 







OR; 



(42) 



This equation is valid in the limit of small particle radius comparing to R. As one can see 
from Eq. fl42]) , optical tweezers measurements can produce a deviation from a real interaction 
potential. 



VII. DISCUSSION 

As shown in this paper, for hard spheres (no mean-field interaction excepting an infinite 
repulsion on contact) in a fluid there are two types of hydrodynamic fluctuation forces: 
(i) the van der Waals forces, mediated by fluctuations of sound waves (U v dw), which is 
similar to the conventional (electromagnetic) van der Waals interaction, and (ii) forces due 
to a thermal drift of particles to the region with a bigger effective mass (/). Despite the 
van der Waals interaction @ formally coincides with the electromagnetic van der Waals 



formula at e — > oo [pq] , existence of such effect in hydrodynamics is not trivial, since a 
calculation on base of the hydrodynamic stress tensor leads to an incorrect conclusion of 
absence of the effect (Section III). The interaction / is also of a fluctuation origin. It has a 
simple mechanical explanation based on that a particle spends more time at a region with 
a bigger mass, since there its velocity is smaller. The proposed hydrodynamic fluctuation 
interactions provide a long range (micron scale) attraction of the order of kT in contrast to 
the conventional van der Waals attraction, which is negligible at the micron scale. 

Charged colloidal particles in an electrolyte, generally, cannot be considered as hard 
spheres, nevertheless, the Coulomb effects modify U v< w weakly, as shown in Section II. The 
interaction I is not modified by the Coulomb effects in a bulk fluid, but can be influenced by 
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an electrostatic restriction of fluctuations in the direction perpendicular to charged confining 
plates. 

The total interaction energy of two particles in a fluid is a sum of the Coulomb repulsive 
part (DLVO) [T|fJ] and the attractive potentials U v dw an d I, plotted in Fig. ^|. A character 
of the resulting interaction depends on particle charges and the Debye length in a very 
delicate way. Formally, the power law fluctuation attractions always win the exponential 
repulsion at some distance Rq. Nevertheless, when Ad is not sufficiently small, the resulting 
potential minimum is far and can be indistinguishable in experiments. Refs. [^Iip2| reported 



an absence of attraction. A lucky choice of electrostatic parameters to observe fluctuation 
hydrodynamic forces corresponds to a close position of the minimum (Rq — 2a) < 2a. 

The hydrodynamic fluctuation interactions is unavoidable and therefore non-interacting 
spheres, strongly speaking, cannot exist in a fluid. The term "hard spheres" means an ab- 
sence of a mean-field (non-fluctuation) interaction excepting an infinite repulsion on contact. 
An interaction between two hard spheres in a dense system of them p3|-p8|, remaining of 
the order of kT, is expected to differ from interaction between two isolated ones, consid- 
ered here, since the hydrodynamic fluctuation forces, due to their non-perturbative nature, 
cannot be reduced to a pairwise interaction. 

The experimentally observed attraction satisfies the fceT-universality condition, i.e., it 
is of the order of k^T in various experiments. This universality is not a trivial property, 
since electrostatic and hydrodynamic mean-field energies are of a few orders of magnitude 
bigger comparing to the room temperature. The /cT-universality enables a strong selection 
of possible attraction mechanisms. Perhaps, a unique candidacy, satisfying this universality, 
is an attraction, mediated by thermal fluctuations of some physical quantities. Since the 
both mechanism proposed above lead to attraction of the order of ZcbT, they are likely 
relevant in interpretation of experimental data on attraction. The next step is to study 
U v dw and / at all R for spherical particles in a bulk fluid and for a confined geometry. A 
contribution of fluctuating surface waves to formation of attraction of particles on liquid-air 
f| and liquid-liquid [|33[ interfaces should increase the attraction effect. 
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FIGURES 
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FIG. 1. The linear chain of particles. Each horizontal segment between particles behaves like an 
elastic harmonic spring. The configuration (a) represents the "natural" positions of two attached 
particles (open circles). The configuration (b) is obtained from (a) by a compression motion of 
attached particles without a destruction of harmonic bonds, (c) is another "natural" position 
obtained from (a) by a destruction of harmonic bonds. 
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FIG. 2. The fluctuation interactions of two cylindrical particles with the axis of the diameter 
length directed perpendicular to two parallel plates, between which they are confined. The in- 
ter-plate distance is of the particle diameter. The dashed line is I(R) and the solid line is the total 
fluctuation hydrodynamic interaction U v jw(R) + 
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